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The upper atmosphere is, as one might expect, an extremely complex system sensi- 
tive to the diurnal changes in the temperature, seasonal and latitudinal varia- 
tions, sunspot activity, and the effects of the solar winds on the geomagnetic 
field. Add to this the mathematical difficulties in solving the gas diffusion 
equilibrium equations for each constituent and the result is a modeling night- 
mare. Nevertheless, from the observations of satellite orbital decay, mass 
spectrometry data, and extreme ultraviolet (EUV) absorption data, Jacchia 
(refs. 1, 2, 3, and M) has developed a series of increasingly accurate models 
which are a careful blend of empirical and theoretical formulae. 


In the most recent Jacchia model, the exospheric temperature is assumed to be 
a function of: 


1) the average and daily variations in the solar flux 

2) the average and three hourly variations in the geomagnetic index 

3) the angle between the position vector and the axis of the unsymmetric 
atmospheric bulge 

4) the angle between the position vector and the geomagnetic pole. 

The exospheric temperature is related to the density by the solution of the dif- 
fusion equilibrium equations for the different constituents of the atmosphere 
as a function of altitude. Other variations are modeled directly as changes in 
the density. They are: 

1) changes due to the semiannual effect 

2) changes due to the seasonal-latitudinal effect. 

The causes for these variations are not exactly known but may be modeled suffi- 
ciently by empirical formulae. The Jacchia model is assumed to be valid over 
the altitude range of 90 to 2500 km. The residuals between the observed den- 
sity from satellite drag observations and the computed densities show the mean 
relative error to be generally less than 10 percent with occasional peak errors 
near 50 percent. 

Although the model recovers most of the important characteristics of the upper 
atmosphere, it is at a great computational expense. The inefficiency of the 
model cannot be reconciled by simply neglecting some of the effects. For exam- 
ple, the 62 standard atmosphere (ref. 5), which assumes that density is solely 
a function of altitude, is grossly inaccurate for altitudes above 150 km (ref. 

6). The major drawback of Jacchia ’s model is the tabular form chosen to repre- 
sent the results of numerically integrating the diffusion equations. The stor- 
age required to implement the table is prohibitive. Walker (ref. 7) has avert- 
ed the storage cost by developing rather complex analytical solutions to the dif' 
fusion equations over a limited altitude range (125-700 km). The analytic ex- 
pressions require several standard function evaluations which hamper the 
computational speed, 

Lineberry, in an effort to reduce the computation time, has assumed that the 
log of the density may be expressed as a truncated Laurent series in tempera- 
ture and altitude. The atmosphere is layered into several altitude bands and 
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the series coefficients of each band are found by pointwise fit to Jacchia's 
tabular results. This layered model together with a more efficient method of 
computing the temperature reduces the computation time five-fold over the 
Walker model and yet maintains the model accuracy and reasonable storage costs. 
In addition, the Lineberry layered expressions are valid over the entire domain 
of the Jacchia model (90-2500 km) and are of such general form that differences 
in the Jacchia models may be easily accommodated. 

This model, which we will call the Jacchia/Lineberry (J/L) model, has already 
been implemented into several orbit propagation software packages (refs. 8 and 
9) but has never been properly documented. The purpose of this report is to 
show the development and concisely define the model, document the accuracy, and 
caution the user of its limitations. 
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The J/L model is most closely fashioned to the Jacchia 71 (ref. 3) model. It 

consists of three basic steps. 

1) Computation of the exospheric temperature requiring the position vector of 
interest, the position of the Sun, the solar flux (F-|q y), and the geomag- 
netic index (Kp) . 

2) Computation of the temperature contribution of the log of the density - 
AZiVt^ requiring the geodetic altitude and temperature. 

3) Computation of the contributions due to the semiannual effect - AUnpc’/^ and 
the seasonal-latitudinal effects - Ajlnog^ requiring the time of the'' year, 
geodetic altitude, and latitude. 

The density, p , is then found from the sum of these contributions 


Znp = (2.1) 


Lineberry assumes as in Mueller (ref. 10) that the temperature contribution has 
the following form 


+ b^z’ * b /z’ (2.2) 

where the coefficients { bj} are selected according to the altitude band of 
interest. The base altitude z' is assumed to be of similar form 


z’ = a^ + a^z + a^/z . (2.3) 

where z is the geodetic altitude. The coefficients { a^} are related to the 
exospheric temperature T^ by the' same function 


^i ^ ^1i ^2i*^“ ^ • (2.4) 

Again, the coefficients selected according to the particular alti- 

tude band of interest. 

The contribution due to the semiannual effect is assumed to be of the form 


^^rip 3 A = + c^z + c^/z)g(t) (2.5) 
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where the coefficients I Cj} are selected according to the altitude band and 
g(t) is a periodic function of time specified in Jacchia’s empirical model. 
Lastly, the seasonal-latitudinal contribution is also assumed to have the famil- 
iar form 


d^/z)p(t)f((j) ) (2.6) 

with band dependent coefficients { dj^} and time dependent function p(t) and 
latitude dependent function f((j) ) specified by Jacchia’s empirical formulae. 

The coefficients and { b j} are considered constants over each altitude 

band and have been determined by pointwise fit of the equations (2. 2-2. 4) to 
Jacchia's tables, being careful to maintain continuity over the altitude band 
boundaries. The width of each band is chosen to maintain three-digit accuracy 
of the assumed expressions and is, therefore, not uniform. The coefficients 
1 Cjj and I djl are also considered constants over each band and are deter- 
mined by pointwise fit of the altitude-dependent parts of (2.5-2. 6) to the 

corresponding altitude-dependent parts of the empirical formulae defined by 
Jacchia . 

An important exception to this general outline is the manner in which hydrogen 
contributions of the density are included. In fitting the coefficients i -j} 
to the Jacchia tables, the contributions due to hydrogen have been removed to 
smooth out the anomalies in the contour plots of the density log versus alti- 
tude and exospheric temperature. The hydrogen contributions may be analytical- 
ly computed directly from the hydrogen diffusion equation by assuming that the 
temperature is a constant. Since the effects due to hydrogen are only signifi- 
cant for altitudes greater than 500 km, where the temperature is almost radi- 
ally isothermal, the assumption is a valid one. 


2.1 INTERPOLATION OF THE JACCHIA TABLES AND EMPIRICAL FORMULAE 

As we have seen, the J/L model uses extensively the interpolation function of 
the form 


f(x) 


= a^ + a^x + a^/x 


(2.7) 


If the left-hand side f(x) is given in some form, such as in Jacchia's tables, 
or by empirical formulae, then the coefficients { a } may be found by inver- 
sion of (2.7). Given the value of the function f(x^) = f^ at three different 
points { x^} i = 1,2,3, the coefficients are then found to be 




(x. 


Ob - ~ 

x^)fl/(x^x^) - 1/(x^x^)) 


( 2 . 8 ) 
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^2 ~ ^^2 “ ^ a^/Cx^x^) 


(2.9) 


= f . - a x^ - a^/x 
1 2 1 3 1 


( 2 . 10 ) 


Equations (2.8-2.10) are to be called the "inverse" relations of (2.7) and will 
be used on many occasions in the following development. 


2-1.1 The Base Altitude Interpolation of Density 

The density at a point is dependent on the altitude and the exospheric tempera- 
ture at that location. If one fixes the temperature at an arbitrary value, say 
T'oo = 600 K, then the density becomes a function of the altitude only. This al 
titude is what we call the base altitude z’ . Equation (2.2) is the interpol- 
ant of the functional dependence of density to the base altitude. The coeffi- 
cients I bjl may be readily determined by the inversion relations along with 
Jacchia’s tables evaluated at the base exospheric temperature T'oo = 600*^K. As 
stated earlier, the effects of hydrogen have been removed from the tables to 
give a smoother interpolant and thus reduce the number of altitude band layers 
in the model. The boundary altitudes of each band were selected in determining 
the coefficients to maintain continuity across the boundaries. The third alti- 
tude was selected so as to give a good fit. The width of each altitude band 
was selected so that the maximum interpolant error was generally less than a 
few percent. The altitude bands and the values of the coefficients { b^} are 
tabulated in table I and table II for the Jacchia 71 and 70 models, respec- 
tively. 


5 
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TABLE I.- JACCHIA 71 VERTICAL PROFILE COEFFICIENTS 


Altitude Band 



•^2 


b 

3 


90-100 

-6.6067 


-1.6401 X 

lO"”* 

1.6968 X 

10^ 

100-110 

-2.2977 X 

10 

-8.2066 X 

10’^ 

9.8734 X 

10^ 

110-140 

-5.4733 X 

10 

6.1437 X 

10~^ 

2.7441 X 

10^ 

140-180 

-3.7147 X 

10 

4.3206 X 

10"^ 

1.4777 X 

iq3 

180-420 

-2.8878 X 

10 

-2.2129 X 

10”^ 

7.2035 X 

10^ 

420-500 

-3.3449 X 

10 

-1.5975 X 

10"2 

1.5545 X 

TO 3 

500-700 

-5.5713 X 

10 

7.7782 X 

io"3 

6.7480 X 

10^ 

700-1500 

-3.8578 X 

10 

-4.8687 X 

10"3 

9.5081 X 

io2 

1500-2500 

-4.1433 X 

10 

-3.8731 X 

10"^ 

2.9930 X 

10^ 
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TABLE II.- JACCHIA 70 VERTICAL PROFILE COEFFICIENTS 


Altitude Band 

^1 



^2 



^3 



90-100 

-2.26064 



-1 .87247 

X 

10"'' 

-3.325619 

X 

lo'' 

100-110 

-2.467081 

X 

10^ 

-7.517851 

X 

10"^ 

1 .087119 

X 

10^ 

o 

1 

o 

-5.856595 

X 

lo"* 

7.759401 

X 

10"^ 

2.967037 

X 

10^ 

140-180 

-3.381609 

X 

lO** 

-9.501784 

X 

10"3 

1 .209134 

X 

10^ 

180-420 

-2.977882 

X 

lo"* 

-2.103046 

X 

10"^ 

8.559544 

X 

10^ 

420-500 

-3.496874 

X 

lo”* 

-1 .404274 

X 

10"^ 

1 .803085 

X 

io3 

500-700 

-5.376797 

X 

10^ 

6.48995 

X 

10"^ 

6.069527 

X 

10^ 

700-1500 

-3.839121 

X 

lo"* 

-4.928746 

X 

10"3 

9.00959 

X 

10^ 

1500-2500 

-4.214804 

X 

10^ 

-3.607654 

X 

10"^ 

3.571183 

X 

10^ 
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2.1.2 The Exospheric Temperature Interpolation of the Base Altitude 

Now if the exospheric temperature is different from that of the base tempera- 
ture T’oo = 600 °K, the base altitude used in (2.2) will be different from that 
of the true altitude. Suppose that at an altitude z and the exosnheric tem- 
perature Too , the density is given by P = P(z,Too). The base altitude is then 
that particular altitude z’ in which the density has the same numerical value 
but at the base temperature T'oo. In other words, z' is defined implicitly by 
the equation 


P(z' , T’^) = p(z, Too) (2.11) 

Equations (2.3) and (2.4) are the explicit interpolants of the above implicit 
equation. The values of the coefficients { a^^} may be determined by selecting 

a particular exospheric temperature tIo and three different altitudes i zj! 

3 - 1,2,3. By equation 2.11 and Jacchia’s modified tables give us three different 
base altitudes {z’} which are implicit functions of the selected temperature 
and altitude. With the values of z’ and z one can invert (2.3) to give the 
coefficients {aj(TOT)} which are implicit functions of the selected exospheric 
temperature. But (2.4) is the explicit interpolant form of these implicit 
functions. By evaluating the coefficients {ajCloo)} for three different 
values of the exospheric temperature 1 looji k = 1,2,3, the inverse of (2.4) 
will give the desired coefficients i . The altitude bands and the value 
of these coefficients are found in table III and table IV for Jacchia 71 and 
Jacchia 70, respectively. 

The interpolant log density versus the exospheric temperature for several alti- 
tudes has been plotted against the Jacchia 71 tables in figure 1. The interpo- 
lant also includes the hydrogen contributions developed in section 2.2. The 
agreement appears to be quite good except for the high altitude, low temperature 
regime in which the hydrogen contribution is significant. The discrepancy can 
be traced to differences in the boundary conditions for hydrogen number density 
in the Jacchia 71 and Jacchia 77 models. The correction in the J/L model to 
include hydrogen assumes boundary conditions as stated in the Jacchia 77 model. 

A comparison to the Jacchia 77 model in figure 2 reflects the similarity in 
the boundary conditions. 
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TABLE III.- JACCHIA 71 BASE ALTITUDE COEFFICIEITTS 


Constant 90-110 km 110-180 km 180-2500 km 


^11 

1.11475 

^12 

1.36100 

^13 

-6.69343 


X 10 3.39245 X 
X 10“^ -5.32690 X 
X 10^ -1.84370 X 


10^ 1.86895 X 10^ 

10”^ 1 .59030 X 10"^ 

10^ -1.17862 X 10^ 


^21 

9.44287 

^22 

7.75000 

^23 

3.31488 

®31 

-5.51954 

^32 

-7.52700 

^33 

3.33882 


X 10“"’ -5.06112 X 

X lO"*^ 2.16963 X 

X 10 8.25561 X 

X 10^ -1.90923 X 

X 10"3 3.23731 

X 10^ 1.02899 X 


10”’' -9.33360 X 10“^ 

10 ”^ 1.34400 X 10“5 

10^ 6.51 163 X 10^ 

10^ -5.47081 X 10^ 

-2.47382 

10 *^ 4.17306 X 10 ^ 
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TABLE IV.- JACCHIA 70 BASE ALTITUDE COEFFICIENTS 



Constant 

qo-110 

km 


1 10-180 

km 


180-2500 

km 

^11 

1.535026 

X 

10^ 

3.86469 

X 

10^ 

1.27264 

X 

102 

^12 

-9.35111 

X 

10“^ 

-7.610145 

X 

10"^ 

4.535789 

X 

10"^ 

^13 

-8.873513 

X 

10^ 

-2.0448485 

X 

10"^ 

-9.268724 

X 

lO** 

a^i 

2.3219^1 

X 

lO"”’ 

-7.287919 

X 

10 "”' 

-3.388665 

X 

10"2 

^22 

4.72682 

X 

10"^ 

3.268459 

X 

10“^ 

-1.339225 

X 

10-5 

^23 

4.43667 

X 

10^ 

9.196106 

X 

10^ 

6.251532 

X 

10^ 

^31 

-7.596 

X 

10 ^ 

-2.158925 

X 

10^ 

4.176991 

X 

103 

^32 

4.58726 

X 

10"”' 

4.417025 



-7.151575 




4.392459 

X 

10^ 

1.136342 

X 

lo"^ 

6.83728 

X 

10^ 

33 
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2 . 1.3 The Semiannual Interpolation 

The variation in the semiannual effect due to changes in the altitude z is 
given in Jacchia (ref. 3) to be 

f(z) = 2.302(5.876 x10’'^z 2*''^31 ^ 0.06328) exp (-2. 868x10“-^ z) (2.12) 

(z in km) 

The factor of 2.302 is needed to convert from the base 10 log used in Jacchia 's 
empirical formulas to the corresponding natural log used in Lineberry’s inter- 
polant ( 2 . 5 ). By evaluating this empirical formulae at three different alti- 
tudes for each altitude band one may determine the coefficients using the in- 
verse relations of the interpolant. The altitude bands and the values of the 
coefficients { Cj^} can be found in table V. The plot of f(z) and its inter- 
polant is shown in figures 3 and 4 for altitudes from 90 to 500 km and 500 to 
2500 km, respectively. The maximum difference between the two amounts to a 
relative error of no more than one percent. 
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TABLE V.- SEMIANNUAL COEFFICIENTS 


Altitude Band C^ C^ 


90-100 

-6.9999 

X 

10-2 

1.4737 

X 

10“^ 

7.8748 



100-110 

-1.2204 

X 

10"^ 

1.1513 

X 

10“^ 

5.3190 



110-140 

-4.6896 

X 

10"^ 

1.3202 

X 

10“3 

7.0920 



140-180 

-1.3067 

X 

10“”' 

1 .6233 

X 

10“^ 

1.2880 

X 

10 

180-420 

-6.5716 

X 

10“^ 

1 .4902 

X 

10“^ 

6.1341 



420-500 

1.0002 



1.5000 

X 

10"*^ 

-2.0940 

X 

10^ 

500-700 

1.6544 



-4.3650 

X 

10“^ 

-3.8535 

X 

10^ 

700-1500 

2.4757 



-1 .0458 

X 

10“^ 

-6.6170 

X 

10^ 

1500-2500 

-8.7290 

X 

10“ 

9.7800 

X 

10“5 

1.788 

X 

10^ 


1M 


Semiannual variation 



J 


Interpol ant 
f(z) 


V'. 

'\ 



'a 



\ 

\ 

\ 

\ 


■\. 


150 


200 

Altitude, km 


250 


300 


350 X 10- 


Figure 4,- Semiannual variation. (500-2500 km.) 
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2.1.4 The Seasonal Latitudinal Interpolation 


The variation in the seasonal-latitudinal effect due to changes in the altitude 
z is given by Jacchia (ref. 3) to be 


g(z) = 2. 302*0. 0l4(z-90)exp( -0.0013(2-90)^'! 
(z in km) (90 < z < l80) 


( 2 . 13a) 


Again, the factor of 2.302 is necessary to convert bases to the corresponding 
Lineberry interpolant (2.6). By evaluating this empirical formula at three dif- 
ferent altitudes for each altitude band, the coefficients of the interpolant 
{d-} may be found by the inverse relations. The altitude bands and coeffi- 
cients are tabulated in table VI. The plot of g(z) and its interpolant is 
shown figure 5 over the valid range of 90 to 180 km. The agreement is quite 
good below 110 km but is less accurate above this altitude because of the width 
of the altitude band. The maximum error results in a relative error in density 
of no more than three percent. An additional band layer could reduce the error. 
But since the error in fitting Jacchia’ s tables is of the same order, the added 
expense may not necessarily increase the total accuracy . 

The Jacchia 70 model (ref. 2) assumes a slightly different form for the seasonal 
latitudinal variation so that 

g(z) = 2.302*0.02(z-90)expf-0.045(z-90)) (2.13b) 

The coefficients fit to this expression are shown in table VII. 


17 



82FM52 



OF POOR 


TABLE VI.- JACCHIA 71 SEASONAL LATITUDINAL COEFFICIENTS 


Altitude Band 


>^2 


D, 




90-100 

8.2812 

-2.8680 X 

10"^ 

1 

o 

o 

X 

10^ 


100-110 

2.4695 X 10 

- 1.1106 X 

lO""* 

-1.3306 

X 

10^ 


o 

o 

5.1205 

-2.4927 X 

10"^ 

-2.1960 

X 

10^ 


lMO-180 

-4.2401 

1 .2570 X 

10“^ 

3.5595 

X 

10^ 
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TABLE VII.- JACCHIA 70 SEASONAL LATITUDINAL COEFFICIENTS 


Altitude Band 

D 

1 


^2 


D 

3 

1 


90-100 

2.4107 X 

10^ 

-1.1142 X 

lO”"* 

-1.2671 X 

10^ 

100-110 

1.5097 X 

10^ 

-6 . 626 X 

10“^ 

-8.1774 X 

10^ 

110-140 

4.3439 


-1.8338 X 

10"^ 

-2.1474 X 

10^ 

140-180 

-1.6246 


3.4375 X 

10”3 

1 .9 4 04 X 

10^ 
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2.2 INCLUSION OF HYDROGEN 

As we mentioned earlier, the contributions due to hydroffen have been removed to 
smooth out anomalies created in Jacchia’s tables. In this section we will 
develop the analytic procedure for includinp; the effects of hydrogen back into 
the J/L model. Since the hydrogen contributions are insignificant at the lower 
altitudes, this development will apply to altitudes of greater than 500 km. 

The barometric diffusion equation governing the number density n(H) 


dn(H) dT 

+ — (1 + a ) + ~dz = 0 (2.14) 

n(H) T R*T 


is where 

T : temperature 

M : molecular weight of hydrogen 

n 

R* : universal gas constant 

g : local acceleration of gravity 

a : thermal diffusion coefficient 
H 

Since at altitudes of z > 500 km the temperature reaches a nearly constant 
value T^, the above equation reduces to 


dn(H) M g 

+ -^dz = 0 (2.15) 

n(H) R 

The acceleration g obeys closely the relation 

g = g (1 + z/R )"^ (2.16) 

o ^ 

where is the acceleration of gravity on the Earth’s surface and R^ is 

the mean Earth radius. Replacing (2.16) into (2.15) and integrating yields 


llnf n(H)l 


+ const. 

R T^(z + Rg) 


(2.17) 
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In ref. 4, Jacchia prescribes the hydrogen number density at an altitude of 500 
km to be 

£nf.n(H)'i^_^QQ = 12.677355 + 66 .544709Tco“^ ( 2 . 18 ) 

which allows us to evaluate the constant of integration. 

The density contribution due to hydrogen is found from 

= M /An(H) 

H n 


or 


jlnp = in(n(H)l ilnfM /a) 
H H ^ 


(2.19) 


with A being Avogadro’s number. 

Combining (2.17) through (2.19) and replacing numerical values for the physical 
constants one arrives at 


= -47.977466 + 66 .544709/Too^ - 7 .006 12x10 ^/Tc» 
+ 7.5572 x10^/(t<»(1 + Z/6378.14)) 


( 2 . 20 ) 
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3.0 THE JACCHIA 71 AND 70 EMPIRICAL MODEL 

The Jacchia 71 (ref. 3) and 70 (ref. 2) model assumes the exospheric tempera- 
ture to be a function of the solar flux, geomagnetic index, and the relative po- 
sition with respect to the diurnal bulge. The exospheric temperature Too used 
in (2.4) has been found to closely obey the following empirical formula. 


T«> = Tl + Tq 

Tq = ATq’Kp + 6TQ*exp(Kp) 

Tl = Ted + R'D) 

= *^C0 ^Vw.7 *^V^10.7 " ^10.7^ 

R = Rq + 6r*Kp 

D = sin”^ + (cos™n - sin"tf )cos'^K/2l 
a = ^Icj) + *^1 
ri = - 4>ol 

T = H + 8 + psin(H + y) 

H - a ~ cLq 

where 

(J) : latitude of given plot 

({>Q : latitude of the Sun 

a : right ascension of given plot 
: right ascension of the Sun 

F- 10.7 : solar flux 10^ Jansky (10“^^ W bandwidth) 

F-]q Y : averaged solar flux 


(3.1) 

(3.2) 


(3.3) 


(3.4) 

(3.5) 

(3.6) 


(3.7) 
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Kp : geomagnetic planetary index 

Kp : averaged geomagnetic planetary index 

and Y, P, B, m, n, R , 5R, , AT , 6T^ , ATp , 6 To 

meters whose values for both models may be round in table VIII. 

The semiannual term g(t) found in (2.5) is given as 


are model para- 


g(t) = 0.02835 + 0 . 3817(1 + 0.467sin($ + 4.l4)lsin(2^ + 4.259) 


(3.8) 


with 


<^> = nt + o.i 9 in{(^ + isin(nt + 6 . 035 ) 


(3.9) 


The Jacchia 70 model represents the semiannual variation indirectly as changes 
in the exospheric temperature. This approach is roughly consistent with the di- 
rect method in Jacchia 71. 

The seasonal-latitudinal terms p(t) and f(c|)) found in (2.6) are for both the 
Jacchia 70 and 71 models defined as 


p(t) = sin(rit + 1.72) (3.10) 

f((j)) = I sin (j)|sin (3.11) 


The time t is in days measured from the beginning of the year (January 1) and 
T\ is the mean motion of the Sun in radians per day. 

In addition, Jacchia corrects for seasonal-latitudinal variations of helium. 
This contribution becomes important at high altitudes and high temperatures 
where helium is a dominant constituent . The variation is expressed in terms of 
a change in the log of the number density of helium. Consequently, the number 
density of helium at the desired altitude and temperature must be known to es- 
tablish this variation. But the J/L model interpolates for the density con- 
sidering all the constituents and not simply helium. Therefore, the J/L model 
neglects this correction. For altitudes z < 500 km the error is negligible. 
At high altitudes and temperatures the relative error can be as large as 100 
percent. 


24 





i\3 














*-3 

o» 

> 

o» 

£> 


a» 

3d 






(D 

^-3 




Q 

30 

O 

3 

3 

oo 

Td 


2 

O 

a 

o 

o 

o 











o 

o 













X 












X 

\ 












— fc 

_a 












o 

o 









G 

O 

o 

-tr 

J= 

o 





o. 

CL 

a 

3 

??! 

X 



X 


1 

1 

1 

(D 

(D 

<D 

H- 



C_ 

c_ 






O'? 

Gf5} 

n 

C+ 



0) 

£B 









M 



3 

D 












Cfl 














isr 






















G 





oo 








PJ 

O 

ro 








1 

+ 

t. 

O 

• 

CO 


'-O 

\0 

O 

o 

OO 

ro 

OO 

CTi 


o 

o 

• 

• 

• 

• 

• 

• 

• 

* 



oo 

3* 

LO 

o 

oo 

^3 

o 

o 

oo 

o 

ro 




H* 




-3r 









CU 







































G 





'njO 








P3 

O 

ro 



oo 





1 

+ 

+ 

O 

« 

00 


OO 

oo 

o 

o 

LjJ 

ro 

oo 

o> 

Sr 

O 

o 

• 

• 

• 

• 

• 

• 

• 

• 

-o 


OO 

X 

UJ 

o 

oo 

'wO 

o 

o 

— X 

o 

-Ji 




H- 




ro 


vD 

oo 






P 






O 

-Cr 



















-O 













O 


>-3 

> 

W 

r* 


< 

M 

l-t 

M 


[?1 

X 

o 

in 

T) 

X 

M 


>-3 

2 


X 

» 

> 

1-^ 

a 

30 

P3 


g 

X 

rq 

r 


X 

3> 

30 

s> 


H 

M 

30 

02 


O O 

=?S 


O 

r=^ 


f© '■'S 

C > 

< 13 


00 
rv) 

z: 

01 


82FM52 


4.0 COMPUTATIONAL CONSIDERATIONS 


.' V.i-i if i' i.vJ 



The empirical formulae given by (3.1) through (3.11) require a number of costly 
standard function evaluations which we would like to avoid if at all possible. 
For example, (3.6) expresses the variation of the temperature with respect to 
the diurnal bulge. Assuming that the position of the point for which we desire 
the density is located by Cartesian coordinates, then four standard function 
calls are necessary to compute the right ascension and latitude used in (3.6). 

If the position of the Sun is also given in terms of Cartesian coordinates then 
we need four more calls to standard functions. Add to this the four calls to 
trigonometric functions found in the calculation of D for a total of 12 stan- 
dard function calls. Since Lineberry's layered expressions allow for more effi- 
cient computation of the density, this diurnal term becomes proportionately a 
major time cruncher. The difficulty is that Jacchia expresses the diurnal 
terms in angular coordinates but the input in most general purpose codes is in 
Cartesian coordinates. The problem can be avoided by expressing D complete- 
ly in terms of Cartesian coordinates without the use of the intermediate angular 
variables . 


An important trigonometric identity 


cos^a = (1 + cos2a)/2 

(4.1) 

sin^a = (1 - sin2a)/2 

(4.2) 


enables us to express the terms in D as follows 

sln’^a = [(1 - sin((|) + (I)q))/ 2]^/2 (4.3) 

cos'^q = [(1 + cos (6 - cj)o) )/2 (4.4) 

and 

cos^|t/2| = (d + cost)/2)^/2 (4.5) 

The value of the terms inside the brackets must always be positive and for this 
reason we can dispense with the absolute value signs appearing in the defini- 
tions of the angular variables. 

Let's examine further the trigonometric term in (4.3). By identity 

sin((}) + ([)q) = sin(j) cos(})q - cos({) sincj)Q (4.6) 
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But the trigonometric terms on the right can be expressed directly in terms of 
the Cartesian coordinates of the point of interest and the Sun. 


sin<b = z/r 

1 

3lnd)o = Zq/Po 

(4.7) 

cos4> = 

+ y^/r ; 

cos4>o = V x^o + y^o/^o 

(4.8) 

The expressions in 

(4.4) and (4.5) 

can be treated in a similar manner. 

A minor 


assumption is necessary to place (4.5) into the desired form. By identity 


COST = cos(H + B)cos(psin(H + y)) ~ sin(H + R)sin(psin(H + y)) (4.9) 

Now the trigonometric functions of H+|8 can be easily reduced to functions of 
the Cartesian coordinates but trigonometric evaluation of the term psin(H + y) 
seems unavoidable. Actually, the magnitude of p is very small, (p = 0.1) 
reflecting the fact that the diurnal bulge is almost symmetric. A Taylor se- 
ries expansion can be used with little error so that 


sin(psin(H+y) ) = psin(H + y) + O(p^) (4.10) 

cos(psin(H + y)) = 1 - (psin(H + y))^/2 + O(p^) (4.11) 

whereas before sin(H + y) is reduced by identity to simple functions of the 
Cartesian coordinates. 

sin H = sina cosa^ - cosa sina^ (4.12) 

cos H = cosa cosa^ + sina sina^ (4.13) 

where 

cosOt = x/Vx^ + y^ ; cosa^ = Xq/^x^q + y^o (4.14) 

sin^t = y/^x^ + ; sin^^Q = y^o • (4.15) 

Since Y and 3 are constants, the trigonometric functions of these terms 
need only be computed once and stored in the database. Also, since in most ap- 
plications the density model is repeatedly evaluated with the time of the year 
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changing little, the trigonometric functions of the position of the Sun need 
not be evaluated frequently. The result is that two calls to SQRT are all that 
should usually be required to evaluate the diurnal term. 


Storing the value of complicated functions of constants or near constants for 
later use is applicable to other computations in the model. For instance, p(t) 
and g(t) found in (2.5) and (2.6), due to their long period, may be considered 
constants over a short period like one day. Also, in many applications, the 
density at some future time must be predicted. Since the short period varia- 
tions in the solar flux and geomagnetic index are next to impossible to predict, 
these values are defaulted to the predicted average values. Here again, the av- 
erage solar flux and geomagnetic index have lengthy periods related to the 1 1 
year solar cycle and may be treated as constants. All these considerations 
have been made in the computational algorithm of the Jacchia/Lineberry model as 
programmed on the flight design system (FDS) (ref. 9). 
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5.0 VERIFICATION OF THE JACCHIA/LINEBERRY MODEL 

The Jacchia/Lineberry model has been implemented in a stand-alone program on 
the UNIVAC 1108-8. Also residing on this computer is the Jacchia 70 model with 
the Walker analytic expressions (J/W) (ref. 11) which allows for a rather 
straightforward comparison. For a meaningful test we have chosen to use the 
Lineberry model parameters fit to the Jacchia 70 model. However, the semi- 
annual variation is represented the same as in the Jacchia 71 model. 


5.1 NUMERICAL COMPARISONS 

For the first comparison, the two models are evaluated at 20 equally spaced 
points, each with the same altitude and all lying in the same plane. The plane 
is oriented so that it is inclined to the equator by M5^ and intersects the 
equatorial plane at a right ascension of 45*^. The epoch date is December 22, 
1977, which places the Sun at its greatest inclination. The average and daily 
solar flux and the average and three hourly geomagnetic index are set at the 
rather quiet conditions of F-jQ.y = 125 and Kp = 2.2. 

To consolidate results, only the numerical average of the densities at the 20 
points are displayed for several different altitudes. The altitudes were cho- 
sen so that they lie just above or below a boundary of the layered atmosphere. 
This enables one to determine if the J/L model exhibits any marked discontin- 
uities. The results of the comparison are shown in table IX. 

The J/L model shows only slight discontinuities across the layer boundaries as 
seen from table V. The boundary at 180 km appears to be the largest, with a 
jump discontinuity in the second digit. l\lhere valid, the J/L and J/W are in 
excellent agreement with discrepancies in the third digit. The largest differ- 
ence is again near the 180 km boundary, but even so the relative error is only 
2.5 percent, certainly within the Jacchia model error. 


TABLE IX.- DENSITY MODEL COMPARISON 


Case 
No . 

Altitude 

(km) 

Averaged 

J/L 

Density (kg/m^) 
J/W 

Relative 
Er ror % 

1 

90 

.344E-5 

N/A 


2 

100- 

.524E-6 

N/A 


3 

100+ 

.524E-6 

N/A 


4 

110- 

.967E-7 

N/A 


5 

110+ 

.965E-7 

N/A 


6 

125 

. 134E-7 

.135E-7 

0.7 

7 

140- 

.384E-8 

.390E-8 

1.5 

8 

140+ 

.384E-8 

.390E-8 

1.5 

9 

180- 

.572E-9 

.558E-9 

2. 4 


29 



OF 


82FM52 


TABLE IX.- Concluded 


Case 

No. 

Altitude 

(km) 

Averaged 

J/L 

Density (kg/m^) 
J/W 

Relative 
Error % 

10 

180+ 

.546E-9 

.558E -9 

2. 1 

11 

420- 

.218E-1 1 

.220E-1 1 

0.9 

12 

420+ 

.218E-11 

.220E-11 

0.9 

13 

500- 

.574E-12 

.582E-12 

1.3 

14 

500+ 

.575E-12 

.582E-12 

1.2 

15 

700- 

.336E-13 

.340E-13 

1.2 

16 

700+ 

.336E-13 

N/A 


17 

1500- 

.581E-15 

N/A 


18 

1500+ 

.573E-15 

N/A 


19 

2500 

.650E-16 

N/A 



To demonstrate the agreement between the models at each point, data from test 
cases 6, 12, 14 and 15 have been plotted showing density versus the angular pa- 
rameter a measured in the plane marking the position of each point. The 
plots are shown on figures 6, 7 , 8, and 9, respectively. Be careful to observe 
the scale of each plot. The agreement is very close for each case. The small 
variations in case 6 are due primarily to the seasonal-latitudinal term. Note 
the significant variations in the density at the higher altitudes, due to the di- 
urnal bulge. 

The last experiment is intended to show that the models agree for different ori- 
entations with respect to the diurnal bulge. In test cases 20 and 21 the Sun 
is placed at the vernal equinox. In case 20 the plane is oriented so that it 
coincides with the plane of the maximum diurnal variation. In contrast, for 
case 21 the plane is chosen to be perpendicular to the maximum diurnal bulge 
plane so we should see little variation. In each case, the altitude is 500 km 
with the F-io 7 = Figures 10 and 11 show the density versus 

angular position for each of the points. The agreement between models is with- 
in the Jacchia model error . 
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Figure 9.- Case 15 density comparison. 
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5.2 TIMING COMPARISONS 

To make an accurate estimate of the computation cycle times, both models have 
been clocked over 100 calls to the algorithm. The cycle time then is the time 
required to evaluate the model 100 times divided by 100. In table X the cycle 
times on the ONIVAC 1110 are shown for both models. The times vary somewhat 
depending on altitude but generally fall within the ranges shown. 


TABLE X.- (DEPUTATION TIMES 



Computation Cycle 

Model 

Time (ms) 

J/W 

3, “I - 3.6 

J/L 

0.5 - 0.7 
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6.0 DENSITY MODEL DIFFERENCES 

The J/L model has been fashioned so that it closely follows the Jacchia 71 
model. The J/L and Jacchia 71 differ in two respects: 

1) The J/L model assumes the hydrogen diffusion boundary condition defined in 
the Jacchia 77 model which is different from that used in the 71 model. 

This difference becomes significant only at high altitudes (z > 500 km) and 
low temperatures (%o 1000° K) . 

2) The J/L model does not correct for seasonal-latitudinal helium variations. 
Only at high altitudes (z ^ 500 km) and high temperatures CHo ^ 500° K) 
does the difference become significant. 

The J/L model differs from the Jacchia 70 model primarily in the manner in 
which the semiannual variation is computed. Since the direct form used in the 
J/L and Jacchia 71 model is roughly equivalent to the indirect form in the 
Jacchia 70 model, the quantitative differences are small. 

The primary quantitative difference between the Jacchia 71 and 77 models is due 
to the form of the geomagnetic effect on temperature. In the 77 model the posi- 
tion with respect to the geomagnetic field is an important feature in describ- 
ing geomagnetic variations in temperature. The difference between models is 
negligible for "quiet" conditions but becomes very large during solar storms 
(large Kp). The Jacchia 77 diurnal variation takes a simpler form than in the 
71 model. The 77 model also suggests methods for determining diurnal varia- 
tions in temperature for each of the constituents. The number density for each 
component is computed according to its temperature. Algorithmically, this repre 
sents a large departure from the 71 model, although numerically the differences 
may be only slight. In such a case, the interpolants in the J/L model must be 
for each constituent instead of the total density. 



. 
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The Jacchia/Lineberry density model is an efficient yet accurate method for com- 
putatiop of the upper atmospheric density. Relative errors between the J/L 
model and Jacchia model are generally less than a few percent. Relative 
errors between the models may become large at high altitudes where the density 
is extremely small. For orbit calculations, such an error is not critical ex~ 
cept possibly for long lifetime studies. The J/L model represents a five- fold 
increase in efficiency with comparable accuracy compared to the J/W model (ref. 
11) used on the real-time computer complex (RTCC) . The modal is of such general 
form that differences in Jacchia models can be readily accounted for. 
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